AD-A102  906  AIR  FORCE  GEOPHYSICS  LAB  HANSCOM  AFB  MA  F/6  22/2 

PREDICTION  OF  MATERIALS  CHARGING  IN  MAGNETOSPHEP IC  PLASMAS. (U) 

MAR  81  AG  RUB IN »  M  F  TAUTZr  K  H  BHAVNANl 
UNCLASSIFIED  AFGL-TR-81-0088 _ NL  _ 


906^’  mm 


AFGL-TR-8 1-0088 

ENVIRONMENTAL  RESEARCH  PAPERS.  NO.  734 


) 


£ 


Prediction  of  Materials 

Charging  in  Magnetospheric  Plasmas 


A.  G.  RUBIN 

M.  F.  TAUTZ 

K.  H.  BHAVNANI 

N.  A.  SAFLEKOS 


19  March  1981 


Approved  for  public  release;  distribution  unlimited. 


SPACE  PHYSICS  DIVISION  PROJECT  7661 

AIR  FORCE  GEOPHYSICS  LABORATORY 

HANSCOM  APB,  MASSACHUSETTS  01731 


AIR  FORCE  SYSTEMS  COMMAND,  USAF 


81  8  17  060 


•*“-v 


This  report  has  been  reviewed  by  the  ESD  Information  Office  (OI)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 


This  technical  report  has  been  reviewed  and 
is  approved  for  publication. 


FOR  THE  COMMANDER 


Cmef  Scientist 


Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Technical  Information  Center.  All  others  should  apply  to  the 
National  Technical  Information  Service. 


^  i 


( 


/I/O  / 

A  l ' 6  I. 


-t  f: 


?1  H 


h 

f' 


a,  r 


/ 


y 


Unclassified 


SECURITY  CLASSIFICATION  OF  this  PACE  (Whan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

T  report  NUMBER 


2  GOVT  ACCESSION  NO 


AFGL-TR-81-0088  „  A  i g  0  Q 


K  T4«b<  /art  ShM«N| 


RF  AO  INSTRUCTIONS 
BEFO  ’  COMPLETING  FORM 


3  RECIP.f- 


BREDICTION  OF  MATERIALS  CHARGING  IN 
MAGNETOSPHERE  PLASMAS, 


7  KUlUBH/.j 

A.  G  J Rubin 
M.  F,  /Tautz 
K.  H.;  Bhavnani 


N.  A.  /Saflekos 


5  T  yPt  o  *■  ‘  * 

Scientific 


ATALOC  NUMBER 


i  pn  •/ 

fnu*r  im 


6  PERFORMING  org  report  h 

ERP  No.  734 


8  CONTRACT  OR  GRANT  NuMBERr#) 


RMtNG  ORGANI  Z  ATtON  NAME  AnO  ADDRESS 

Air  Force  Geophysics  Laboratory 
Space  Physics  Division  (PHK) 
Hanscom  AFB,  MA  01731 


10  PROGRAM  element,  p 
AREA  8  WORK 


76610703 

F7TT51F 


MEnT.  PftfWiCT,  TA» 
UNIT  NUWBE  RS 

Cli- 


n  CONTROLLING  OFFICE  name  and  aooress 

Air  Force  Geophysics  Laboratory 
Space  Physics  Division  (PHK) 
Hanscom  AFB,  MA  01731 


1  j 


19  MaraiMfl81 


13.  NUMBER  of  pages 
20 


4  monitoring  ~a£Tn  CY  NAM^  ADORESSf/J^lV/eranr  from  Controlling  Office) 


'5  SECURITY  CLASS.  (of  this  report) 

Unclassified 


<6.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 


‘1 


/  , 


^  / 


L 


>  r  f. 


17.  DISTRIBUTION  STATEMENT  (of  fhe  •b»fr«cr  «nf«r»d  In  Block  20.  if  different  from  Report) 


V 


18  SUPPLEMENTARY  NOTES 


19  KE  Y  WORDS  f  Continue  on  reverse  aide  if  necessary  end  Identify  by  block  number) 

Spacecraft  charging 
Materials 

Geosynchronous  plasma 
Secondary  emission 


o  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

*  Basic  materials  charging  in  magnetospheric  plasmas  may  be  predicted 
based  on  material  properties  that  define  capacitance,  conduction,  and  emission 
for  the  surface,  and  on  plasma  species,  temperatures  and  densities.  This 
report  summarizes  results  obtained  using  MATCHG,  a  program  that  embodies 
the  basic  formulations  of  the  NASCAP  program  without  the  3 -dimensional  (3-D) 
complexities.  Plasma  temperature  and  material  property  variation  effects  on 
surface  potentials  are  considered,  A  graphical  procedure  for  predicting  volt¬ 
age  response  from  the  general  no  bulk  conductivity  characteristics  is  presented 


DO 


1473 


Unclassified 


SECURITY  CLASSIFICATION  of  THIS  PAGE  »'h#n  Dot*  Fn 


woio  ,  v 


."pi  in mm 


A ccesilm 


DiStr • 
/•vr*  1 


»v— t !  <*n/ 


Cc-Ios 


1  ;.r)U/or 


Contents 


1.  INTRODUCTION 

2.  MATERIALS 

3.  CHARGING  AT  LOW  TEMPERATURES 

4.  EFFECT  OF  BULK  CONDUCTIVITY 

5.  GRAPHICAL  AND  LINEARIZED  REPRESENTATION  OF 

CHARGING  EQUATIONS 

6.  CIRCUIT  EQUATIONS 

7.  GRAPHICAL  SOLUTION 

8.  NUMBERS  FOR  T  =  5.  6  CASE 


Illustrations 


1.  MATCHG  Potentials  for  Low  Temperatures 

2.  MATCHG  Potentials  for  Intermediate  Temperatures 

3.  MATCHG  Potentials  for  High  Temperatures 

4.  MATCHG  Potentials  vs  Backing  Voltage  (Bulk 

Conductivity  Included) 

5.  MATCHG  Potentials  (Bulk  Conductivity  Included) 


Illustrations 

(j.  MATCHG  Potentials  —  Teflon  (Bulk  Conductivity  Included)  14 

7.  MATCHG  Circuit  Diagram  17 

8.  Net  Current.  J  vs  V  From  MATCHG  17 

9.  Calculated  V  vs  Time  From  MATCHG  IP 

Tables 

1.  Material  Properties  7 

2.  Double  and  Single  Maxwellian  Plasma  Cases  10 

3.  Single  vs  Double  Maxwellian  Results  11 

4.  MATCHG  Potentials  vs  NASCA  P  Quasisphere  12 

5.  MATCHG  Potentials  vs  SC1-3  14 


4 


Prediction  of  Materials 
Charging  in  Magnetospheric  Plasmas 


I.  INTRODUCTION 

It  is  often  useful  to  obtain  a  quick  approximation  of  the  potential  to  which  a 
material  would  charge  in  a  geosynchronous  orbit  plasma.  If  one  does  not  con¬ 
sider  multidimensional  effects  and  the  effects  of  photoelectrons,  a  computation 
of  the  materials  effects  alone  is  possible. 

This  calculation  is  carried  out  by  means  of  a  computer  code  called  MATCHG 
(Materials  Charging).  Given  a  plasma  environment  of  specified  electron  and 
proton  temperature  and  number  density,  the  code  calculates  the  potential  to  which 
a  material  will  charge  in  this  plasma.  The  principle  of  the  method  is  to  solve  for 
the  potential  for  which  there  is  zero  net  current  to  the  material  surface.  Included 
in  the  current  balance  are  incident  electron  and  proton  (ion)  currents,  secondary 
electrons  from  both  incident  electrons  and  ions,  backscattering,  and  current 
leakage  to  a  substrate  through  bulk  conductivity.  In  the  formulation  which  includes 
bulk  conductivity,  the  thickness  of  the  material  must  be  specified.  The  secondary- 
emission  yield  curves,  as  a  function  of  energy,  are  calculated  in  the  code. 


(Received  for  publication  17  March  1981) 


2.  MATERIALS 


The  materials  for  which  calculations  were  carried  out  are  those  typically 
used  on  satellite  surfaces  (see  Tabic*  1).  "SOLAR"  refers  to  solar  cell  cover 
glass  material;  "WHITEN"  to  a  nonconducting  white  ”  YELOWC"  to  a  ;  •  !1  w 

cotidiit"  . ng  paint;  "BLACKC"  to  a  black  conduct’  ...  paint. 

Hie  MATCHG  program  charges  a  surfac*  by  nu'nt's  .  f  material  !* 

•  i . i: i : . *  i •  >:u.  ;,r-  used  in  NASCA  P.  1  :i  ■  si-*  1.  ,.a<-  *  its  -.Cl;  .•.-Mtnr.*-  ' 

the  equilibrium  potentials  for  a  surface  immersed  in  a  spr  <1  plasma  envir;  n 
ment.  The  full  3-0  complexities  of  more  realistic  models  are  neglected. 

In  this  report  we  present  the  main  results  from  MATCHG.  Graphical  and 

linearized  techniques  are  developed  for  approximating  the  results  of  the  MATCHG 

equations.  Runs  were  made  with  the  plasma  temperature  T  (T.  =  2T  )  as  the 

e  g  i  _  ^  e 

independent  variable.  The  densities  were  fixed  at  =  n.  =  10  M  .  The  poten¬ 
tials  were  computed  in  three  ranges; 

1.  Low  temperatures  (see  Figure  1) 

2.  Intermediate  temperatures  (see  Figure  2) 

3.  High  temperatures  (see  Figure  3). 

Runs  were  made  to  calculate  the  potential  for  a  two-Maxwellian  plasma  case 
and  the  corresponding  single  Maxwellian  plasma  case.  The  parameters  specify 
three  types  of  plasmas:  hot,  cold,  and  average.  These  are  shown  in  Table  2. 

The  MATCHG  results  are  given  in  Table  3  for  Teflon  and  aluminum. 

The  intermediate  temperature  cases  were  done  for  all  of  the  materials  cur¬ 
rently  available  in  MATCHG  (see  Figure  2).  Note  that  materials  BLACKC  and 
YELOWC  are  identical,  that  is,  they  are  represented  in  the  program  by  the  same 
material  parameters.  Also  the  materials  WHITEN  and  Kapton  are  sufficiently 
near  to  YELOWC  and  BLACKC  that  the  final  potentials  are  the  same  (to  three 
decimal  places). 


3.  CHARGING  AT  LOW  TEMPERATURES 

Spacecraft  charging  is  most  significant  for  plasma  temperatures  greater  than 
1  kv.  It  is  of  interest  to  study  charging  at  lower  plasma  temperatures  as  well. 

For  nonsunlit  conditions,  Figure  1  shows  the  equilibrium  potentials  for  a 
number  of  materials,  both  metals  and  dielectrics. 


1.  Katz,  I.  et  al  (1977)  A  Three  Dimensional  Dynamic  Study  of  Electrical  Charg¬ 
ing  in  Materials.  NA§A  CR-135^56, 
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Table  2,  Double  and  Single  Maxwellian  Plasma  Cases 


1) 

Hot  Plasma 

CD 

o 

o 

0> 

c 

n . ,  -  0.  06 

il 

ne2  0. 14 

n>2  °‘  14 

T  ,  -  1200  eV 

el 

T  -  2000  eV 
il 

T  „  =  10000  eV 
e2 

T.41  -  1G000  eV 

i2 

1 -Maxwellian  version: 

n  =  n.  =  0.  20 
e  i 

T  =  8000 
e 

T.  =  13000 

i 

2) 

Cold  Plasma 

n  ,  =  0.  85 
el 

n.  =0.  25 

1 1 

ne2  1  °-25 

n.2  =  0.85 

T  -  =  1000  eV 
e  1 

T.,  =  300  eV 

1 1 

T  0  =  4000  eV 
e2 

T  -  9000  eV 

1 -Maxwellian  versions: 

n  =  n.  =  1.  1 
e  i 

T  =  1800 
e 

T.  =  8000 

i 

3) 

Average  Plasma 

nel  =°-9 

nu  =  °-9 

ne2  *0'6 

n.2  =0.6  | 

T  ,  =  1500  eV 
el 

T.,  =  427  eV 
ll 

Te2  =  4400  eV 

T,2  =  10000  eV 

1 -Maxwellian  versions: 

n  =  n.  =  1.  5 
e  i 

T  =  2800 

T.  =  7000 

e 

i 

=  7000 


V  (VOLTS) 


Te  (keV) 


Figure  3,  MATCITG  Potentials  for  High  Temperatures 

At  extremely  low  temperatures,  up  to  5  eV,  all  materials  charge  negatively 
to  approximately  the  same  potential.  Between  10  and  100  eV,  secondary  emis¬ 
sion  properties  become  important  and  the  potentials  of  different  materials  peak 
between  -20  V  for  Teflon  to  -115  V  for  gold.  Between  about  20  eV  and  200  eV 
plasma  temperature,  different  materials  charge  to  positive  voltages.  The  posi¬ 
tive  voltages  are  all  less  than  5  eV,  most  around  2  eV. 

In  the  temperature  region  between  10  and  several  hundred  eV,  materials 
potentials  are  either  a  few  volts  positive  or  tens  of  volts  negative. 

These  results  are  significant  in  the  design  of  electron  spectrometers  where 
the  potentials  produced  by  low-energy  electrons  impinging  on  dielectric  materials 
is  important. 

4.  KFFKCT  OF  BM.K  CONDI  CTIVITY 

MATC'HG  includes  the  efrect  of  the  bulk  conductivity  current.  In  order  to 
calculate  this  current,  the  potential  difference  across  the  charging  surface  must 
be  known  that  is,  the  voltage  underlying  the  surface  is  needed. 


Table  3.  Single  vs  Double  Maxwellian  Re  sults 


Initial 

Final 
Potent  lal 

Rat  ii  »s 
(MAX  ' 

J  .Ie<  f  t  ‘  ,n 
(  u t rent 

R :  O  :  ■  ■ 

(MAX 

Phis:'; a  Type 

Material 

<k«V) 

DMf.MAXI 

M0 

DPI. MAX  t 

MAX  -  h.  • 

Aiummu:.; 

-0.  00  ^ 

l .  n  i 

-0.  4  7'*  ^ 

DH1. MAN  - 

A lu: :  -m;:  . 

.4  / 

-0.  I'll  J 

MAX  - 

Alunmu:::; 

""  } 

0.  88 

J.o- 

DHL  MAX  -  aw 

Alum  mu::; 

-1.-4  f 

-2,00  > 

MAX  -  4-old 

Aluminum 

-o.i 

1.90 

1  .  10 

DRLMAX  -  cold 

Aluminum 

-o.ohoo  > 

-1.  14  J 

MAX  -  hot 

Teflon 

0.94  ! 

-0.  4  79  | 

1.11 

DRI-MAX  -  hot 

T  eflon 

-2.90  ' 

-0. 420  J 

MAX  -  ave 

Teflon 

1 

-0.  00071  'i 

l.Ofi 

} 

1. 00 

DBLMAX  -  ave 

Teflon 

*0.00007 /  1 

-2.00  J 

MAX  -  cold 

Teflon 

+0. 00127  | 

0.  92 

, 

1.  10 

DBI.MAX  -  mid 

To  flon 

+0.  0012  8  J 

-1.14  J 

To  make  a  comparison  to  a  NASCAP  satellite  model,  the  potential  on  the 
underlying  conductor  would  be  taken  as  the  backing  plate  voltage.  Table  4  shows 
the  shade-side  potentials  for  quasispheres  in  a  fixed  sun  environment  and  the 
corresponding  MATCHG  voltages,  including  the  bulk  conductivity  current.  There 
is  good  agreement  for  all  of  the  materials. 

In  Figure  4  we  plot  the  MATCHG  potentials  vs  the  backing  plate  voltage  for 
the  case  of  Teflon  in  the  nominal  T  =  10,  20  keV  plasma.  This  figure  shows  that 
the  minimum  voltage  |  -v|  occurs  when  the  backing  voltage  is  0. 

In  Figure  5  we  show  the  lower  limit  MATCHG  potentials  (backing  plate  voltage 
set  to  0)  for  a  number  of  different  materials.  The  plasma  temperature, 

T  (T.  =  2T  ),  is  the  independent  variable  and  the  densities  are  fixed 

€  i  ^0 

(ne  =  n.  =  10  M  J.  It  is  evident  that  those  materials  with  large  bulk  conductiv¬ 
ities  (WHITEN,  SI02)  are  depressed  relative  to  the  others,  and  SOLAR,  with  a 
small  conductivity,  charges  up  higher,  that  is,  the  larger  the  conductivity,  the 
easier  it  is  for  charges  to  leak  through  the  material  and  reduce  the  voltages. 


Table  4.  MATCHG  Potentials  vs  NASCAP  Quasisphere 


Material 

NASCAP 

Shade 

Potential 

Quasisphere 

Conductor 

Potential 

MATCHG 

Backing 

Plate 

Final 

Potential 

Teflon 

-8894 

-7125 

-7100 

-8800 

Kapton 

-9210 

-7378 

-7380 

-9280 

SOLA  R 

-7820 

-(>284 

-(>2  GO 

-7820 

SI02 

-358 

-279 

-280 

-3G0 

WHITEN 

-2997 

-2395 

-2390 

-2990 

Plasma:  T  -  10,  T.  =  20  keV,  n  -  n.  =  10(l  M 
e  i  e  i 

Sun  is  in  direction  (1,  1,  1) 


Figure  4.  MATCHG  Potentials  vs  Backing  Voltage  (Bulk  Conductivity 
Included) 


i 


2 


Figure  5,  MATCHG  Potentials  (Bulk  Conductivity 
Included) 


Lower  limit  curves  of  this  type  would  be  directly  applicable  to  cases  where 
the  spacecraft  conductor  does  not  charge  and  the  surface  lies  in  the  shade  (and 
the  surface  conductivity  is  negligible).  For  example,  we  can  compare  the  num¬ 
bers  obtained  for  the  SCATHA  model  for  experiment  SC  1-3  (which  lies  in  the 
shade)  to  the  MATCHG  results,  as  shown  in  Table  5.  The  agreement  is  very 
good. 

When  the  bulk  conductivity  current  is  included  in  the  computations,  the  steady- 
state  potentials  are  determined  by  a  balance  between  the  current  arriving  from 
the  plasma  and  that  escaping  through  the  surface  to  the  backing  plate.  The  final 
voltages  reached  will  thus  depend  on  the  plasma  density  (which  is  related  simply 
to  the  plasma  current).  This  was  not  the  case  before  when  the  bulk  conductivity 
was  neglected  (see  Figures  1  and  2).  In  Figure  6,  we  show  the  large  changes 
that  take  place  in  the  MATCHG  potentials  when  the  plasma  density  is  varied  by 
factors  of  10  and  0.  1  (for  the  case  of  a  Teflon  surface). 

One  can  get  a  feeling  for  the  effect  of  the  ion  temperature  of  the  plasma  on 
the  MATCHG  potentials  from  Figure  6.  The  plot  shows  the  change  in  the  voltage 
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curve  that  occurs  when  we  use  the  relation  =  T.  instead  of  the  nominal 

=  2T\.  There  is  only  a  small  increase  10  percent)  in  the  equilibrium  poten¬ 

tials. 

An  alternate  formulation  for  the  secondary  emission  was  added,  as  a  test 

mode,  to  MATCHG.  In  this  new  version,  the  Sternglass  yield  equation  is  used 
2 

at  low  energies 


9  m  _'r*N  t 

6/6m-xt'-  (x-E,  r.  ' 


.  3 


and  the  Wall  formula  at  higher  energies 
*=kx-°-725  . 


By  imposing  continuity  of  slope  and  value  for  the  above  two  equations  we  set 
the  transition  energy  x  =  3  and  the  ratio  K/6m  =  1.  54.  The  MATCHG  potentials, 
from  runs  using  this  test  mode,  are  shown  in  Figure  6.  The  test  equations  result 
in  a  slight  increase  in  the  number  of  secondaries  and  thus  the  equilibrium  poten¬ 
tials  are  decreased,  as  indicated  in  the  figure . 


5.  GRAPHICAL  AND  LINEARIZED  REPRESENTATION  OF 

CHANGING  EQUATIONS 

We  solve  the  circuit  equations  representing  the  charging  calculation  under 
the  assumption  that  the  source  of  current  is  a  linear  function  of  the  source  volt¬ 
ages.  This  procedure  gives  a  fairly  good  approximation  to  the  MATCHG  results  — 
the  general  features  of  the  charging  process  are  obtained  (including  final  equilib¬ 
rium  voltage  and  time  constant,  which  depends  on  conductivity,  thickness  of 
material,  etc. ). 

6.  CIRCUIT  EQUATIONS 

The  charging  calculations  that  are  done  in  the  MATCHG  program  can  be  rep¬ 
resented  by  the  circuit  shown  in  Figure  7.  The  equation  describing  the  voltage 
on  the  capacitor  is,  from  the  diagram, 

2.  Dekker,  A.J.  (1958)  Solid  State  Physics,  Academic  Press,  New  York  (i:279. 

3.  Wall,  J.A.  et  al  (1977)  Results  of  Literature  Search  on  Dielectric  Properties 

and  Electron  Interaction  Phenomena  Related  to  Spacecraft  Charging, 
Proceedings  of  the  Spacecraft  Charging  Technology  Conference.  AFGL- 
nC-77 -005L  AD  A045T59. 
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dV 

dt 


i-(V-VB)oA/d  x 
_  = 


Ke  A/d 

o 


(dj  -  (V 


VB)a) 


(1) 


where  j  -  i/A  is  the  new  current  density  from  the  source  (it  does  not  include  the 
bulk  conductivity  current).  MATC'HG  calculates  j,  which  is  a  complicated  func¬ 
tion  of  the  plasma  environment  and  the  properties  of  the  surface  material.  In 

Figure  8,  we  show-  the  j  j(V)  dependence  for  Teflon  in  a  plasma  specified  by 

-  3 

T  5  keV,  T.  -  G,  n  0.  1  cm  ,  n .  =  0.  15.  It  is  evident  that  the  i  vs  V  curve 

e  *  •  i  •  e  1  J 

is  almost  linear, 


e  ’i 

Thus  it  will  be  a  fairly  good  approximation  to  take 


j=.i0  +  J'V  . 


where  J'  is  some  average  slope.  For  simplicity,  we  could  let 


(2) 


J'  =  - 


J,. 


V(o  =  0) 


(3) 


where  j  =  initial  net  current  density,  V(cr  =  0)  =  equilibrium  potential  when  a  =  0 
which  gives  the  straight  line  indicated  in  Figure  8  (V(ct  =  0)  =  0  and  j  <  0  so  that 
J’  <  0).  Substituting  Eq.  (2)  into  Eq.  (1)  gives 


d?  =  i<7  (djo  +  vB)  + v(k7“<j'  "°/d))  =  A  +  BV  . 


(4) 


where  we  define 


A  =  Kl-(d3o  +  VBCT) 

O 


B  =  -  a/d)  <  0  . 

o 


Equation  (4)  is  easily  solved.  We  write  it  as 


(5) 


(6) 


dV  _  _ 

(A  +  BV)  "  dt 


Integrating  both  sides  of  Eq.  (7)  gives 
4  In  (A  +  B  V)  =  t  +  C  . 


(7) 


(8) 
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k  ^DIELECTRIC  CONSTANT 
€0*8.85xIO"12 
d  =  THICKNESS  OF  MATERIAL 
A  =  AREA  OF  SURFACE 
a  *BULK  CONDUCTIVITY 


dV.l c_ 
dt  C 

lC=  '  -  'R 


EQUATIONS 
FOR  CAPACITOR 

charging 


Figure  7.  MATCHG  Circuit  Diagram 


V  (VOLTS) 


Figure  8.  Net  Current  J  vs  V  From  MATCHG 
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The  integration  constant  C  can  be  evaluated  by  letting  t  =  0  in  Eq.  (8),  that  is, 


r 

i 

l 


C  =  ^-ln  (A  +  BV  )  ,  (9) 

±3  O 

where  Vq  =  V(t  =  0).  Putting  Eq.  (9)  into  Eq.  (8)  and  solving  for  V  leads  to 


In  Figure  9  we  plot  V  from  Eq.  (10)  (for  the  case  Vq  =  0)  vs  the  time  and  the 
corresponding  results  obtained  from  MATCHG.  The  predicted  curve  is  slightly 
higher.  The  time  constant  for  charging  up  can  be  taken  as 


r 


B 


Ufo 

d 


1 

(J'  -  a/d) 


Keo 


1 

(J'd/a  -  1) 


(11) 


Check:  when  J*  =  0  (constant  source  current),  we  have 


which  is  what  we  would  expect  from  elementary  circuit  theory.  The  final 
equilibrium  potentials  can  be  obtained  by  setting  =  0  in  Eq.  (4),  or  from 
Eq.  (10)  by  setting  t  =  oo,  since  the  exponential  term  goes  to  0  (B  <  0). 


ve  ■  v<->  -  -  -f 


(J'  -  a/d) 


or  using  Eq.  (3), 


=  V(a  =  0) 


Check:  when  a  =  0  we  recover  the  result  that  V  =  V(a  =  0). 

e 


(12) 


Figure  9.  Calculated  V  vs  Time  From  MATCHG 


7,  GRAPHICAL  SOLUTION 


As  time  increases,  the  source  current,  i,  decreases  and  the  conductivity 
current,  i^,  increases  in  proportion  to  the  voltage  buildup  on  the  capacitor.  At 
equilibrium,  the  current  to  the  capacitor,  i^,,  has  decreased  to  0  and  thus 
i  =  i^,  or  in  terms  of  current  densities  we  have 


j=jo  +  J’Ve=JR  *  ?<Vo-VB> 


Solving  the  condition  for  Ve  again  yields  the  Eq.  (12).  One  could  solve  the 
above  equations  graphically  by  finding  the  intersection  of  the  straight  lines  rep¬ 
resenting  j  and  jp  on  the  (j  ,  V)  plot  (see  Figure  8).  In  fact,  an  exact  solution 
would  be  determined  by  the  intersection  of  the  line  and  the  MATCHG  j(V)  curve. 
The  effect  of  changing  a,  d,  and  is  clear  from  this  construction.  One  just 
changes  the  slope  and  intercept  of  the  line  and  finds  the  new  intersection  point 
giving  Ve. 
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